o 

(N 

u 

< 



O 
U 

I 

O 



> 

^' 

O 
en 



X 



Stochastic Gravitational Wave Background generated by Cosmic String Networks: 
Velocity-Dependent One-Scale model versus Scale-Invariant Evolution 
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We compute the power spectrum of the stochastic gravitational wave background generated by 
cosmic string networks described by the Velocity-Dependent One-Scale (VOS) model, for a wide 
range of macroscopic and microscopic parameters. The VOS model — which has been shown to 
provide an accurate macroscopic description of the evolution of cosmic string networks — is used to 
demonstrate that cosmic string networks are unable to rapidly attain scale-invariant evolution after 
the transition between the radiation and matter eras. However, in computations of the stochas- 
tic gravitational wave background, it is often assumed that the networks experience scale-invariant 
evolution throughout cosmological history. We demonstrate that this assumption leads to an un- 
derestimation of the amplitude and broadness of the peak of the spectrum, that may consequently 
lead to inaccurate observational constraints on the cosmic string tension. 



I. INTRODUCTION 



The production of cosmic string networks as remnants 
of symmetry breaking phase transitions in the early uni- 
verse is predicted in a wide variety of grand unified sce- 
narios. Initial attempts to model cosmic string networks 
[l|-y| were based on the assumption that a single length- 
scale is sufficient to describe their dynamics. These one- 
scale models were later ameliorated by treating the aver- 
age root-mean-square velocity of the network as a dynam- 
ical variable, while maintaining the one-scale assump- 
tion. This Velocity dependent One-Scale (VOS) model 
[J, [5| provides a quantitative description of the string 
network throughout its evolution and it has been thor- 
oughly tested and calibrated using numerical simulations 
of cosmic string network dynamics [6-9] . For a general- 
ization of the VOS model for networks of domain walls 
or p-brane of arbitrary dimensionality, see |10l - ll5l | . 

Cosmic string interactions are crucial to the dynam- 
ics of cosmic string networks. These interactions re- 
sult in the formation of closed cosmic string loops that 
detach from the long string network. Cosmic string 
loops oscillate under the effect of their tension and de- 
cay by emitting gravitational waves. This emission oc- 
curs throughout the cosmological history and gives rise 
to a Stochastic Gravitational wave background (SGWB) 
with a broad frequency range |16l - [l9| . Different regions of 
the stochastic gravitational wave background power spec- 
trum may be probed using current and upcoming astro- 
physical experiments: direct gravitational wave detectors 
— eith er g round-based (Advanced LIGO [20], Advanced 
Virgo [111 and KAGRA jH]) or space-borne (evolved 
LISA/NGO [Hi and DECIGO [HDrmilsar timing exper- 
iments (Parkes [25] and European f26'| pulsar timing ar- 
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rays and NANOGRAV [2^); small-scale fluctuations and 
B-mode polarization of Cosmic Microwave Background 
[Hi (Planck ^ and CMBpol ^); and big bang nu- 
cleosynthesis [31|. There is thus the prospect either for 
the detection of the SGWB generated by cosmic string 
networks or for the tightening of current constraints on 
the cosmic string tension. 

The shape and amplitude of the stochastic gravita- 
tional wave background power spectrum is highly de- 
pendent on both the large-scale properties of the net- 
work and the size and emission spectrum of the cos- 
mic string loops. Previous estimations and computa- 
tions of the SGWB spectrum [32iJ41i| generated by cos- 
mic string networks were performed using the one-scale 
model, thus implicitly assuming that the networks expe- 
rience scale invariant evolution throughout cosmological 
history, and that the scaling parameters suffer a sudden 
change at the radiation-matter transition. This, how- 
ever, does not give an accurate picture of network evo- 
lution: during the transition between the radiation- and 
matter-dominated eras, the network adapts slowly to the 
changes in the underlying background dynamics. In this 
paper, we study the SGWB generated by cosmic string 
networks described by the VOS model. We will find that 
the assumption of scale-invariant evolution during the 
radiation- and matter-dominated eras has indeed a sig- 
nificant impact on the shape and amplitude of the spec- 
tra and may, thus, result in inaccurate constraints on the 
cosmic string tension. 

This paper is organized as follows. In Sec. |TT1 we review 
the VOS model for cosmic string network dynamics, and 
provide a brief description of network evolution through- 
out cosmological history. In Sec. IIIIl we describe the 
emission of gravitational waves by cosmic string loops, 
and present a method for computing the SGWB power 
spectrum that results from the superimposition of the in- 
dividual emissions by many cosmic string loops. In Sec. 
llVi we compute the stochastic gravitational wave back- 
ground generated by cosmic string networks for a wide 



variety of parameters and analyze the effect of consid- 
ering a realistic radiation-matter transition on the shape 
and amplitude of the spectrum. We then conclude in Sec. 
lYl 



II. MICROSCOPIC COSMIC STRING 
NETWORK EVOLUTION 

The world-history of an infinitely thin and featureless 
cosmic string may be described by a two-dimensional 
worldsheet, 



obeying the Nambu-Goto action 

S = —ji 



CpUy 



(1) 



(2) 



where a = 0, 1, u^ {u^) are timelike (spacelike) co- 
ordinates that parameterize the cosmic string world- 
sheet, and /i is the cosmic string tension. Here, g = 
det(go-iy), 9(71^ = 9avx"a,x\ is the induced metric, and 
gai, is the metric tensor of the background. In a flat 
3-|-l-dimensional homogeneous and isotropic Friedmann- 
Robertson- Walker (FRW) background, the line element 
is given by 



Ss — gavdx" dx^ = a^{rf) {d^rj — dx • dx) 



(3) 



where a represents the cosmological scale factor, t and 
r] = J dt/a are respectively the physical and conformal 
times, and x is a 3- vector whose components are comov- 
ing cartesian coordinates. 

By varying the action in Eq. ([2]) with respect to x"' , 
and imposing temporal-transverse gauge conditions. 



V, 



0, 



the string equations of motion can be written as [42 



where 



e + 2Hyi^e = 0, 

^/2 \ 1/2 



,'\ ' 



1-x 



C2 



(4) 

(5) 
(6) 

(7) 



H = d/a = Ha, H is the Hubble parameter, and dots and 
primes represent partial derivatives with respect to the 
conformal time and to the spacelike parameter u = u^, 
respectively. 



A. Velocity-Dependent One-Scale Model 

The cosmic string equations of motion (Eqs. (O and 
^) can be averaged in order to obtain evolution equa- 
tions which describe statistically the large-scale evolu- 
tion of cosmic string networks. This might be done by 



assuming that the network is roughly homogeneous on 
sufficiently large scales — and thus it may be described 
by an unique lengthscale — and by treating the root- 
mean-square (RMS) velocity of the network as a dy- 
namical variable. This model — the Velocity-Dependent 
One-scale (VOS) Model [J, Q — provides a quantitative 
description of the large-scale evolution of cosmic string 
networks from early to late cosmological times (see Refs. 
[ij, [l^ H^l for a more general approach) . 

The total energy in cosmic strings, E, and the RMS 
velocity of the string network, v, are defined respectively 
as 



(8) 



E = fia{rj) / edu 



(x^) 



Jx^edu 
J edu 



(9) 



The average string energy density, p, is proportional to 
Ea~^. By using Eqs. ([5]) and ([SUHl) and neglecting, for 
the moment, the energy loss caused by loop production, 
one obtains 



dp 
'dt 



2Hp{l + v^) = 0. 



(10) 



Cosmic string interactions are a key ingredient in the 
evolution of cosmic string networks. When two cosmic 
strings collide, they may exchange partners and inter- 
commute. This process has important consequences for 
the dynamics of a cosmic string network, since it leads to 
the generation of cosmic string loops which detach from 
the long-string network. These loops have a finite lifes- 
pan, leading to an energy loss by the string network. A 
cosmic string network then has mainly two constituents: 
long strings that stretch over large scales (larger than the 
horizon), and smaller cosmic string loops. The large-scale 
properties of the long string network may be described by 
an unique scale — the characteristic length L — defined 
as 



_ M 



(11) 



where p is the long string energy density. The rate of 
energy loss caused by loop production may be written as 



a 



dp 



loops 



~-p 

cv— . 
L 



(12) 



where c is a phenomenological parameter which char- 
acterizes the efficiency of this energy loss mechanism 
(this term needs to be added to Eq. (|TU)) ). Numerical 
simulations indicate that, for ordinary cosmic strings, 
c = 0.23 ± 0.04 is a good fit both in the matter and 
radiation eras [5]. 

Using Eqs. ([TU)l - (112p . one may find the evolution equa- 
tion for the characteristic length, L, 



2— 
dt 



2H+—]L + dv. 



(13) 



Here, we have introduced the damping lengthscale, 



2H 






(14) 



which accounts, not only for the damping caused by the 
background expansion, but also for the effect of frictional 
forces caused by the interaction of cosmic strings with 
other cosmological components (encoded in the friction 
lengthscale if). Throughout this paper, we shall assume 
that If = oo, so that i~^ = 2H. 

The evolution equation for v might be obtained by dif- 
ferentiating Eq.(ini) and using Eqs. ([SEl) 



dv 
'di 



(1 - v^) 



(15) 



Here fc is a dimensionless curvature parameter that de- 

gends, in general, on v (for a detailed definition of k see 
|[lj|). In 0, the following ansatz was suggested 



k{v) = ^ (1 



«2) (l + 2V2f;''') - 



1-8-0^ 
"8^ 



(16) 



Note that the energy loss due to loop production is in- 
cluded implicitly in Eq. (J15p through the dependency of 
the curvature term on L. 

In order to also take into account the effect of the gravi- 
tational back-reaction on the dynamics of the long string 
network, one needs to add an extra term to the right- 
hand side of Eq. p^ : 



dL 
'dt 



= ATGuv'^ 



(17) 



where F ~ 65 [a, [la 01 is a parameter which describes 
the efficiency of the emission of gravitational waves by 
long strings. 



B. The Scale-Invariant Regime 



is an attractor solution of the VOS equations (Eqs P^ 
I15[)l. in the case of decelerating power-law expansion of 
the universe — with a oc t^ and < /3 < 1. (For a 
more detailed discussion of the scaling solutions of cosmic 
string networks, see |14l[l5ll4g|.) 

This scaling solution is only attainable for a constant 
expansion exponent /3 and, therefore, such a regime can 
only be maintained deep into the radiation or matter 
eras. Moreover, the values of the scaling constants are 
different in the radiation and matter epochs. One would 
then expect the network to experience two different scale- 
invariant regimes separated by a transitional period (trig- 
gered by the radiation-matter transition) during which it 
is not in linear scaling. This issue has been discussed in 
detail in [49] , and the authors found that this transition is 
indeed long-lasting, especially for small values of c. This 
is illustrated in Fig. [1] where the cosmological evolution 
of ^ and V is plotted as a function of the scale factor for 
a cosmic string network in a 3 + 1-dimensional FRW uni- 
verse, for different values of c. Notice that, due to the 
recent acceleration of the expansion of the Universe, the 
matter epoch might not be long enough for the network 
to attain a linear scaling regime before being diluted away 
by the accelerated expansion. 



III. STOCHASTIC GRAVITATIONAL WAVE 
BACKGROUND 

Cosmic string loops are copiously produced through- 
out the cosmological evolution of cosmic string networks. 
Once these loops detach from the long string network, 
they oscillate relativistically under the effect of their ten- 
sion, and decay by emitting gravitational waves. As a 
result, cosmic string networks are expected to give rise 
to a stochastic gravitational wave background spanning 
a wide range of frequencies and having a characteristic 
shape [l6l.ll7|. 



It has been demonstrated numericallyjTL l45l l46p and 
using a variety of semi-analytical models 0,[3j|43,|4^ that 
cosmic string networks may evolve towards a linear scal- 
ing regime during which the characteristic length grows 
proportionally to the particle horizon. In such condi- 
tions, the string energy density remains a constant frac- 
tion of the background energy density in both matter- 
and radiation-dominated epochs. Cosmic strings net- 
works are thus not expected — unlike domain walls and 
magnetic monopoles — to dominate the energy density 
of the universe at late times. 

A linear scaling regime of the form 



L = ^t 



and 



with 



e = 



k{k^ 



4/3(1 - /3) ' 



constant . 



I k 1-/3 
k + c P 



(18) 



(19) 



A. Cosmic String Loop Emission 

In this paper, we shall assume that loops are born with 
a size which is a fixed fraction of the characteristic length 
of the string network at that time (t^) 



lb = aL{tb) , 



(20) 



where I — E/fi is the physical length of the cosmic string 
loops, E is the energy of the loops, a is a constant param- 
eter and the subscript b refers to the instant of cosmic 
time when the loops were produced. Realistically, one 
does not expect all the loops produced at a given time to 
have precisely the same length. Instead, the distribution 
of the sizes of the loops formed at the time it, is expected 
to have a peak around lb- If the width of the distribu- 
tion of loop sizes is not very large, then assuming that all 
the loops have the same size at the moment of formation 
should be a good approximation. 
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FIG. 1: Evolution of ^ (top panel) and of the RMS velocity v 
(lower panel) of a network of cosmic strings with G/i — 10^^ 
in a fiat 3 + 1-dimensional FRW universe, as a function of 
the scale factor a, for various values of c (see label). The 
cosmological parameters were set to JIao = 0.728, Q,roh^ — 
2.47 X 10"^ and h = 0.704, which correspond to the WMAP 7- 
year data combined with the baryon acoustic oscillation data 
from the Sloan Digital Sky Survey and determinations of Ho 
using the Hubble Space Telescope [501 . The scale factor at the 
time of radiation and matter equality is a^q, while its value 
at the present time was set to unity. 



Loops emit gravitational waves in a discrete set of fre- 
quencies, which are mainly determined by the length of 
the loop 



f -^ 



(21) 



where j is the harmonic mode and fj is the corresponding 
frequency. The energy emitted in gravitational waves by 
a cosmic string loop in each of the modes is 



where 



dt -^'^^ 



m— 1 



(22) 



(23) 



r ~ 65 [la . |44| is a constant that characterizes the effi- 
ciency of the emission mechanism, and q is the spectral 
index. The value of q depends on the shape of the loops: 
it should be q « 2 for kinky loops, while we should have 
q « 4/3 for cuspy loops [51]. There is another effect 
tliat needs to be taken into consideration when consid- 
ering loop emissions: gravitational back-reaction. It has 
been demonstrated, using both field theory simulations 
and analytical studies [5^ |53| that the effect of grav- 
itational back-reaction is to damp the higher frequency 
modes while leaving the lower modes virtually unaffected. 
A cut-off may, then, be introduced to the summation in 
Eq. (I^5|) . so that emission modes with j > Ug are not 
considered. In this case, Eq. (P^ remains valid with oo 
replaced by an integer cut-off n^. In Sec. lIVBi we shall 
study the dependence of the SGWB spectrum on Us and 
determine the values that this parameter should take. 
Cosmic string loops, then, loose energy at a rate 



f--"^ 



(24) 



and, as a consequence, their length decreases. As a mat- 
ter of fact, we have that 



l{t)=aL{tb)-TG^l{t^tb), 



(25) 



for t > tb, assuming, for simplicity, that the loops do 
not break up into smaller loops due to self-intersection. 
Loops created in this manner would have a faster decay 
rate than that of their progenitors and they would emit 
gravitational waves with larger frequencies. 



B. Spectral Density of Gravitational Radiation 

The amplitude of the stochastic gravitational wave 
background is often quantified by the energy density in 
gravitational waves, pcW) per logarithmic frequency in- 
terval in units of critical density (/Ocrit), 



f^GW — 



1 dpGw 

Pcrit dlogf 



(26) 



In order to find the spectral energy density of gravita- 
tional waves at a given time i, one needs to take into con- 
sideration all the contributions from gravitational waves 
that have a frequency / at the instant t, emitted by cos- 
mic string loops created between the time of formation of 
the network, tf, and t. Taking into account the redshift 
of th e g ravitational wave frequency (/ oc a^^), one has 
that fsil 



df 



(t) 



27r 



*.t'^(*' 



a 



it) 



X 



X / dln(l,t')lh(2TTfl^ 
Jo \ a{t') 



(27) 
(28) 



where n{l,t')dl is the number density of cosmic string 
loops with lengths between I and I -\- dl at time t, and 
h is a function that describes the spectrum of radiation 
emitted by a loop. If one assumes a discrete spectrum 
(as in Eqs. (|22p and (|23p l. this function may be written 
as (13 



n{lj{t'),t') has contributions from all preexisting loops 
that have physical lengths lj{t') at time t'. Determining 
the times of creation (i|,) of the loops that contribute to 
a given frequency at any given time t' is essential to cor- 
rectly computing n(lj{t'),t'). Given these instants, one 
has that 






(29) 



where z = {a{t) / {a{t')) {2tt ft) . Therefore, one finds that, 
at the present time. 



f^Gw(/) 



PcritOOo/ Jtf 



ta 



' dta^t')Y,jP,nil,it'),t') , 



j=i 



(30) 

where the subscript is used to refer to the value of the 
corresponding parameter at the present time. Here we 
have also defined 



iUt'),t')^j2''Uti),ti) 



a{t') 



(33) 



where we have taken into account the fact that the loops 
are diluted away by the background expansion. The num- 
ber density of loops created at the instants t\ may be 
computed using Eq. 



n(;,(i^),i^) = 



We then have that 



duc 
~dt 



dt 
ll 



(34) 



W) 



2j ajt') 
f 0,0 



(31) nil,{t'),t') 



1 



vitl) 



which is the physical length that the cosmic string loops 
should have at each instant t' in order to emit, in the har- 
monic mode j, gravitational waves that have a frequency 
/ at the present time. Given this dependence of I on i', 
the relevant components of the loop distribution func- 
tion, that may contribute to r^GW^ at a given frequency 
f,a.Teniljit'),t'). 



C. Number density of Loops 

In order to compute ficwj it is necessary to correctly 
compute the loop distribution function nil jit'), t'). Let 
nc{t) be the total number density of loops that have been 
formed as the result of intercommutation between the 
time of formation of the string network and time t. If 
one starts by evolving the VOS equations numerically 
— in order to compute the characteristic length, L, and 
the RMS velocity, v, for a discrete set of cosmic times 
between tf and to — the rate of loop production (per 
comoving volume) may easily be determined: 



duc 



c V 



(32) 



This expression may simply be obtained by dividing the 
total energy density that is lost by the string network due 
to loop formation (Eq. (fT2)l ) by the energy of each loop 
at the moment of creation (recall that, for simplicity, all 
loops produced at a given time are assumed to be created 
with the same physical length and thus the same energy). 
Note that, after formation, loop size shrinks as a con- 
sequence of gravitational radiation emission. Therefore, 



dL I 




. (-^U^,.+rGf,aL^{tl)\ait') 

(35) 
Note that, in performing this calculation — contrary to 
what is generally done in the literature ^37i . i3S-i4l| — 
we did not assume the network to be in a linear scal- 
ing regime. As we have discussed in Sec. IIIBl once 
the radiation-matter transition is triggered, for a real- 
istic expansion history and realistic values of c, the cos- 
mic string network is not expected to remain in a scale- 
invariant regime. As we shall see in the rest of this paper, 
considering a realistic radiation-matter transition has a 
significant impact on the shape of the spectrum at low 
frequencies. 



IV. THE EFFECT OF THE 

RADIATION-MATTER TRANSITION ON THE 

SGWB SPECTRUM 

The stochastic gravitational wave background pro- 
duced by a cosmic string network is, at any given fre- 
quency, the result of the superimposition of the contri- 
butions of all the loops that generate, throughout the 
cosmic history, gravitational waves with that frequency 
at the present time. The loops created during the radi- 
ation epoch are the dominant contributors to the high 
frequency range and they produce a flat spectrum. On 
the other hand, larger loops, created during the matter 
era, generate a peak in the low frequency portion of the 
spectrum. 

Although the general shape of the SGWB spectrum 
generated by a cosmic string network is maintained for a 
wide range of parameters, its amplitude and the charac- 
teristics of the peak may vary greatly. In this section, we 



shall compute the SGWB spectrum generated by cosmic 
string networks for a wide range of parameters. We per- 
form these computations using two different models for 
the evolution of cosmic string networks: the VOS model 
described in Sec. Ill Al (Model I) and a model whose un- 
derlying assumption is that cosmic string networks expe- 
rience scale-invariant evolution throughout cosmological 
history (Model II). In Model II, the required linear scal- 
ing behavior is achieved by determining the values of the 
scaling parameters during the radiation and matter eras 
(using Eq. (|19p l. and by assuming that v and ^ change 
abruptly between these values, in a step-like manner, at 
the time of radiation and matter equality. Note that the 
approximation made in Model II is quite common in the 
calculation of the SGWB spectrum generated by cosmic 
string networks. Computing the SGWB power spectra 
for both these models will allow us, not only to assess 
the differences in the shape of the spectra that result 
from the variation of cosmic string parameters, but also 
to study the impact of an inappropriate modeling of the 
radiation-matter transition on the shape of the spectra. 
Throughout this paper, we assume a 3 -I- 1-dimensional 
FRW background containing matter, radiation and a cos- 
mological constant. The cosmological parameters take 
the values I^ao = 0.728, flrO^J^ = 2.47 x IQ-^ and 
h = 0.704 [50|, while F = 65. The fiducial value of the 
energy loss parameter is c = 0.23. 



A. Dependency on a and Gfj, 

The shape and intensity of the SGWB spectrum is af- 
fected by the value of the loop size parameter a. There 
is presently no consensus in the literature as to what 
value should a take. Several studies suggest that a 
should be smaller than the gravitational back-reaction 
scale a < TGfj,, while others suggest scales much closer 
(albeit 1 to 3 orders of magnitude smaller) to the char- 
acteristic lengthscale of the network (see [5J-[6J|). There 
are also studies which indicate that string loops might 
be formed with lengths similar to the string thickness 
|65l - l67l |. and others suggesting that the loops generated 
at any given time may be associated with two or more 
fundamentally different scales |68l.l69j. Given the uncer- 
tainties associated to the value of the loop size parameter, 
we shall compute the SGWB power spectrum for a wide 
range of values of a. 

There are two different regimes to take into consid- 
eration when analyzing the effects of a on the shape 
and intensity of the spectrum. Large loops, with sizes 
above the gravitational back-reaction scale (a > FG/x), 
have a lifetime larger than the Hubble timescale {H~^). 
Consequently, in this regime, cosmic string loops per- 
sist for a significant amount of time before disappearing 
completely, emitting gravitational radiation throughout 
their lifetime. One may then expect the density of grav- 
itational waves to be larger in the case of loops with a 
larger loop size parameter a. This effect is shown in 



Fig [51 where the SGWB power spectrum generated us- 
ing Model I is represented for different values of a. Fig. 
[2] also shows the analogous SGWB power spectra (char- 
acterized by the same Gjjl and a) in the case of Model 
II. The analysis of the plots shows that Model II pre- 
dicts lower and narrower peaks, when compared to Model 
I. This happens because, for realistic values of c, once 
the radiation-matter transition is triggered, the param- 
eter ^ suffers a slow increase, while v decreases slowly 
(as is evident in Fig. [Ij. Assuming that this variation 
occurs abruptly introduces a discontinuity in n{l{t'),t') 
and leads to an underestimation of the number of loops 
chopped from the long-string network during the matter- 
dominated era. Additionally, assuming an instantaneous 
change in the scaling constants causes an overestimation 
of the size of the loops at the moment of creation during 
the matter epoch and the consequent overestimation of 
the lifetime of the loops. Although this has the effect of 
increasing the peak of the spectra, it is not sufhcient, in 
general, to counteract the other effect described above. 
At low frequency, there is also a slight difference in the 
spectra that is caused by the recent acceleration of ex- 
pansion caused by dark energy. The effect of dark energy 
is as one might expect: it causes a decrease of the ampli- 
tude of the spectrum at low frequencies, because of the 
decline in the number of loops produced that is caused 
by the accelerated expansion (during such a phase L ex a 
and u — > {14|). 

In the small loop regime (a < FG/x), cosmic string 
loops have a lifespan smaller than the cosmological time 
scale and may, therefore, be assumed to radiate their 
energy immediately after being chopped. It is straight- 
forward to realize that, in this limit, the shape of the 
spectrum is independent of a. However, in this case, 
the effect of changing a is to change the frequency of 
the emitted gravitational waves. The spectrum, then, 
merely suffers a shift towards higher frequencies as a in- 
creases, as is clearly illustrated in Fig. [31 As to the 
differences between the spectra generated using Models 
I and II, they are mainly caused by the overestimation 
of ^ during the matter epoch. As matter of fact, a loop 
created with a size I radiates gravitational energy with 
frequencies / > 2/1. Loops created during the matter 
era by a network that remains artificially in a linear scal- 
ing regime during the transitional period would be larger, 
and, as a consequence, the SGWB spectrum generated by 
these networks exhibits a shift towards lower frequencies. 

The effects of varying the string tension, Gfi, on the 
amplitude of the spectra may be explained similarly. 
First of all, one may once again analyze the large loop 
and small loop regimes separately. The spectral density 
depends on Gfi directly (see Eq. ([50)) ). but also indi- 
rectly since the lifetime of the loops depends on G^. In 
the small loop regime, the effect of varying Gfi on the life- 
time of the loops is negligible — since loops radiate their 
energy effectively immediately after formation. The am- 
plitude of the spectrum, then, varies proportionally to 
string tension (Q,QT^h^ ex G^). If the loops are large. 
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FIG. 2: The SGWB spectrum, SIgw'i , as a function of the 
frequency, /, for various values of a > VGji,. The sohd lines 
represent the spectra created by a cosmic string network de- 
scribed by the VOS model, while the dashed lines represent 
the spectra obtained using Model II. Here, we assumed that 
G\ji = 10~^ and c = 0.23, and only the fundamental mode of 
emission was considered (ns = 1). 



FIG. 4: The SGWB spectrum, JIg-w^^, as a function of the 
frequency, /, for various values of G[i. The solid lines rep- 
resent the spectra generated using Model I, and the dashed 
lines represent the spectra obtained using Model II. The loop 
size parameter was set to a = lO"'^, the loop-chopping effi- 
ciency was set to c = 0.23, and only the fundamental mode 
of emission was considered (n^ = 1). 
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FIG. 3: The SGWB spectrum, Qcwh^, as a function of the 
frequency, /, for various values of a < FG^. The solid lines 
represent the spectra generated using Model I, while dashed 
lines represent the spectra obtained using Model II. As in Fig. 
[21 we have assumed that Gfj, = 10~^ and c = 0.23, and ris — 1. 



ho-wever, the effect of decreasing G/i is to increase the 
lifetime of the loop. This causes a less pronounced de- 
crease of the SGWB spectrum amplitude and a shift of 
the peak to-wards higher frequencies. These effects are 
clearly illustrated in Fig. |4l -where the SGWB spectra 
generated using Models I and II are represented for var- 
ious values of G^, and for a — 10~^. Note also that 
the differences bet-ween Models I and II persist for dif- 
ferent values of Gfi and they are quite significant in the 
mid-to-small frequency range. 



B. Dependency on n^ and q 

The spectrum of gravitational radiation emitted by 
cosmic string loops depends both on the cut-off on the 
number of frequency modes, Us, and the spectral index, q. 
As we shall see in this section, the shape of the stochastic 
gravitational -wave spectrum is highly dependent on the 
values of these parameters. 

Figs. [5] and [6] sho-w the spectral density of radiation 
generated by a cosmic string net-work -with a large loop 
size parameter {a = 10~^) for different values of Ug and 
q. The inclusion of higher order modes of emission causes 
a shift of the peak to-wards higher frequencies. This shift 
is, in both cases, accompanied by a decrease of the max- 
imum amplitude of the spectrum and by a broadening of 
the peak. These effects may easily be inferred from Eqs. 
(22) and ([^5)1 : they are simply a consequence of -weighing 
in higher order frequency modes. One might also deduce 
from these expressions that the spectrum generated by 
cosmic string net-works -with cuspy loops {q = 4/3) is sig- 
nificantly more affected than that of string net-works "with 
kinky loops {q = 2). In both cases, these differences are 
more significant for larger values of rig. However, one 
may notice, by analyzing Figs. [5] and [6l that the spec- 
tra seem to converge to a fixed shape for large values of 
Us ■ Although in the case of kinky loops it is sufficient to 
consider modes up to rig — 10^, for cuspy loops the cut- 
off should be some-what larger (n^ = 10'') (these results 
seem to be in agreement -with those found in Ref. |41]). 

When one considers a small loop size parameter, the 
picture is some-what changed. One may see in Figs. [7]and 
|5] that the effect of including higher order modes is qual- 
itatively the same. Ho-wever, the shape of the spectrum 
is significantly affected for cuspy loops {q = 4/3). In 
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FIG. 5: The SGWB spectrum, Jlcwft'^, as a function of the 
frequency, /, for different values of n^. The sohd hnes rep- 
resent the spectra generated by a cosmic string network de- 
scribed by VOS model (Model I), while the dashed lines are 
spectra obtained using Model II. Here, we have considered 
kinky loops (with q = 2) with a loop size parameter a — 10~^, 
and we have assumed that G^ = 10~^ and c = 0.23. 
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FIG. 7: The SGWB spectrum, JIgw^-^, as a function of the 
frequency, /, for different values of n^. The solid lines rep- 
resent the spectra generated by a cosmic string network de- 
scribed by the VOS model (Model I), while the dashed lines 
are obtained using Model II. Here, we considered kinky loops 
(with q = 2) with a loop size parameter a = 10~*, and we 
have assumed that Gjj, = 10~^ and c = 0.23. 
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FIG. 6: The SGWB spectrum, SIgw^ , as a function of the 
frequency, /, for different values of n^. The solid lines rep- 
resent the spectra generated using Model I, while the dashed 
lines are the spectra obtained using Model II. Here, we con- 
sidered cuspy loops (with q — 4/3) with a = 10~^, and we 
have assumed that Gfi — 10~^ and c — 0.23. 



FIG. 8: The SGWB spectrum, Qawh^, as a function of the 
frequency, /, for different values of n^. The solid lines rep- 
resent the spectra generated using Model I, while the dashed 
lines are obtained using Model II. Here, we considered cuspy 
loops (with q — 4/3) whose size is determined by a = 10~*, 
and we have assumed that Gfi — 10~^ and c — 0.23. 



this case, the broadening of the peak is very significant 
and the spectrum develops a hump as it decreases to- 
wards the constant portion. Moreover, for a < FG/i, one 
needs to consider values of n^ that are about one order 
of magnitude higher to achieve convergence of the shape 
of the spectrum. Notice also that the spectra exhibits 
small "bumps" in the low frequency range. These are 
a consequence of the fact that loops emit gravitational 
waves in a discrete set of frequencies. Recall that the 
gravitational waves emitted by the j-th harmonic have a 
minimum frequency of / = jfo (where /o is the frequency 
of gravitational waves emitted at the present time in the 



fundamental mode) and, thus, this mode does not con- 
tribute to the spectral amplitude of gravitational radia- 
tion for / < j/q. Each successive harmonic comes into 
play at / — jfo, thus causing a sudden change in the 
slope of the SGWB spectrum. Note that these disconti- 
nuities are more evident for a small loop size parameter 
because, in this case, the values / = jfo coincide with the 
peak region (whose slope is quickly varying). Note also 
that these "bumps" would not be present if one would 
consider the more realistic case of the loop size following 
a distribution of values around aL{tb) (instead of consid- 
ering only one loop size). 



By analyzing Figs. [5] to [51 one may conclude that the 
differences between the spectra generated using Models 
I and II are still present for large values of Ug and for 
q = 2 (kinky loops) and q = 4/3 (loops with cusps). 
For large values of a, the differences are essentially of 
the same magnitude as in the case of n^ = 1. For small 
a, however, the maximum amplitude of the spectrum 
generated using Model I surpasses that obtained using 
Model II for the Ug > 1 case (despite the fact that the 
amplitude of the peaks was very similar for n^ = 1). 
This happens because, in the case of Model II, the peak 
is situated at very low frequencies, close to /o- This is 
the region that suffers a more accentuated decrease as a 
result of considering larger values of n^ . 



C. Dependency on c 

String-theory-motivated scenarios suggest that funda- 
mental strings (F-strings) and 1-dimensional Dirichlet 
branes (D-strings) may grow to macroscopic scales and 
play the cosmological role of ordinary cosmic strings (see 
e.g. [204221 )■ These cosmic superstrings — unlike ordi- 
nary cosmic strings — have a reconnection probability, P, 
that may be significantly smaller than unity: for instance, 
for D-string intersections 0.1 < P < 1, while for collisions 
between F-strings the expectation is 10"'^ ^ ^ ^ 1 [23- 
Hence, in D-D or F-F collisions, there is a non- vanishing 
probability of the strings merely passing through each 
other without interaction. 

It is straightforward to realize that, as a result of the 
reduction of the reconnection probability, the amount of 
loops created throughout the evolution of the network 
is also significantly reduced. The energy-loss mechanism 
of these networks is, thus, less efficient, and the loop 
chopping parameter is expected to scale as c ex P. For a 
fixed value of Gfj,, weakly interacting cosmic superstring 
networks have, as a consequence, a larger string energy 
density and thus an enhanced SGWB amplitude. 

As a matter of fact, during the radiation era, weakly 
interacting string networks with c <^ 1 attain a linear 
scaling regime characterized by |49|, 



^ = V2d 



V2 



(36) 



with S = (7r/12)5. Using Eqs. pi^ and ([351), we 
may then conclude that, for weakly interacting networks, 
n{l,t) ex P~^ and / ex P~^, during the radiation era. As 
a result, the amplitude of the flat portion of the SGWB 
spectrum should scales as 



n 



GW 



cxP" 



(37) 



As to the matter era, the fact that the network is not 
experiencing scale invariant evolution makes quantitative 
analysis more difficult. Fig. [T] shows that the effect of 
varying c on u and ^ is non trivial. On the one hand. 



for c ^ 1, the scaling law ^ ex P is approximately main- 
tained. However, for larger values of c, the increase in 
^ is somewhat smaller and, as consequence, the decrease 
on the number of loops created that occurs as a result of 
the radiation-matter transition is less accentuated. Con- 
sequently — although for very weakly interacting net- 
works (P < 1) Eq. ([37|) is stiU valid — for P - 1 the 
decrease of the amplitude of the peak of the SGWB spec- 
trum is less accentuated. 



V. CONCLUSIONS 

A cosmic string network is not expected to maintain a 
linear scaling regime throughout the cosmic history. Dur- 
ing the radiation-matter transition, the network needs to 
adapt to the changing background conditions, entering 
a phase in which the parameters ^ and v are not con- 
stant. Moreover, the matter epoch is not long enough for 
the network to be able to attain scale-invariant evolu- 
tion before the universe becomes effectively dark-energy- 
dominated. Therefore, it is not realistic to assume that 
cosmic string networks remain in a scaling regime af- 
ter the transition between the radiation- and matter- 
dominated eras. 

Most of the computations of the SGWB power spec- 
trum available in the literature assume that the net- 
work experiences linear scaling evolution, with the scaling 
parameters changing abruptly at the time of radiation- 
matter equality. In this work, we have assessed this as- 
sumption by studying in detail the spectra generated by 
cosmic string networks that undergo a realistic cosmolog- 
ical evolution. This was done for a wide range of parame- 
ters and the results were compared with those of models 
in which the strings are assumed to remain in a linear 
scaling regime. We found that the linear scaling assump- 
tion has significant impact on the shape and amplitude of 
the spectrum, in particular, in the low-to-mid frequency 
range. This range corresponds to the peak of the spec- 
trum which has, in general, an amplitude that is a few 
orders of magnitude larger than that of its flat portion, 
and may then be within the reach of future experiments. 

In Fig. [HI some examples of SGWB spectra generated 
using Model I that are within the sensitivity range of the 
European Pulsar Timing Array (EFTA) experiments |74| 
are plotted. In all these examples, the spectra obtained 
using Model II would not seem to be detectable. This 
clearly illustrates the importance of accurately modeling 
cosmic string network dynamics, when computing the 
Stochastic gravitational wave background. The power 
spectra computed under the assumption that the cosmic 
string networks remain in a linear scaling regime have 
narrower peaks, characterized by a lower amplitude and 
a larger slope, for a wide range of parameters. This im- 
plies that simplified models of string evolution may lead 
to inaccuracies in the computation of the constraints on 
cosmic string tension Gfi using gravitational wave de- 
tection experiments. In particular, for those experiments 
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FIG. 9: Some examples of SGWB spectra generated by cosmic 
string networks (see label for details) that would not seem to 
be within the detection range of the 5-years (pink area) and 
10-years (gra y a rea) European Pulsar Timing Array (EPTA) 
experiments [TJI, when scale invariance is assumed (dashed 
lines). In these examples, the equivalent spectra obtained 
using the VOS model (solid lines) are within the detention 
range of the EPTA experiments. 



that have a sensitivity window in the low frequency range, 
the constraints obtained using the linear scaling assump- 
tion are weaker than those obtained using the VOS model 
for the evolution of the cosmic string network. 
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